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ABSTRACT 


The viscous, radiating hypersonic flow past an axisymmetric 
blunt body is analyzed based on the Navier-Stokes equations, plus a 
radiative equation of transfer derived from the Milne -Eddington dif- 
ferential approximation. The fluid is assumed to be a perfect gas with 
constant specific heats, a constant Prandtl number, P, of order unity, 
a viscosity coefficient varying as a power, u), of the temperature, and 
an absorption coefficient varying as the first power of the density 
and as a power, n, of the temperature. 

The gray gas assumption is invoked, thereby making the absorp- 
tion coefficient independent of the spectral frequency . Limiting 
forms of the solutions are studied as the freestream Mach number, 

M, freestream Reynolds number, R, and the temperature ratio across 
the shock wave, T = eM , go to infinity, and as the Bouguer number, 

Bu, and the density ratio across the shock wave, e = (y-l)/(v+l), go 
to zero. 

The method of matched asymptotic expansions is used in the anal- 
ysis, and it is shown that there is a far-field precursor, composed of 
two regions, in which the fluid mechanics can be neglected for all 
practical purposes but are included for completeness. Also found is 
a near-field precursor in which the temperature rises to its shock 
structure order of magnitude. The high temperature in this region, 
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due to radiation transfer, causes the next region, the outer region 
of the shock structure, to be thicker by OCr^) than the outer region 
•without radiation. Next follows the middle and inrier regions of the 
shock structure, with solutions differing only in the constants of 
integration from those found for the non-radiative problem. A radia- 
tion relaxation region between the shock structure and shock layer is 
now needed in order that radiative equilibrium can be obtained. The 

shock layer is unchanged from the non-radiative problem, except for 

1 ). 

the addition of I = T , and has two or three regions, depending on the 
ratios of the viscous to inviscid terms. 
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Chapter I 
INTRODUCTION 

The purpose of this paper is to isolate the important parameters 
of a viscous, radiating hypersonic flow past a blunt body of arbitrary 
shape. Much confusion and discussion have existed as to the importance 
of the radiative term's contribution to the energy equation in the 
various flow regions from the precursor field immediately ahead of 
the shock wave to the body surface, i.e., references [13] and [17J. 

In order to isolate the parameters, and the role they play, and 
to make the problem tractable, the gas is assumed to be perfect with 
an absorption coefficient that is independent of spectral frequency 
(the gray gas assumption). The gray gas assunption is a serious one 
in that spectral detail in the problem is completely omitted, which 
could allow the gas to range from a high absorption at one frequency 
to a low absorption at another frequency, thereby having absorption 
properties ranging from optically thick to optically thin in the same 
flow region. On the other hand, the gray gas assumption yields the 
problem analytically tractable, thereby allowing zhe functional be- 
haviors of the various flow quantities from the freestream to the body 
to be found. If spectral detail were included, provided it were 
known, only numerical results could be hoped for, which would give no 
indication as to how the various flow quantities interact with each 
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other, Heaslet and Baldwin [lj . Figure 1 gives the results of 
Heaslet and Baldwin for the strong shock case. The parameter k/k' 
corresponds to the inverse of the Boltzmann number that is used in the 
present analysis . 

Most analyses up to now have assumed that the various regions of 
flow were either optically thick (the Rosseland approximation) or op- 
tically thin (the Planck approximation), Emanuel [V], Jische J12J 
Burggraff [2], and Olfe [15] . The thick and thin approximations are 
the lower and upper limits, respectively, of the photon mean free path 
and therefore are even more restrictive than the gray gas assumption. 
These two approximations have the additional serious drawback that one 
of the two conditions for similarity is eliminated. More specifically, 
the thick and. thin approximations identically satisfy the 
transfer equation, thereby eliminating the Bouguer number, inverse 
characteristic optical length, and its importance as a similarity 
parameter . Although the Bouguer number appears in the energy equation 
it appears only in the ratio of Bouguer number to Boltzmann number. 

On the other hand, the condition for "full similarity" is the equiva- 
lence of the ratio Bu/Bo and the Bouguer number. In the present treat- 
ment, the full coupling of the fluid mechanics and radiative transfer 
are used, with the details of the analysis being the determining factor 
as to whether or not the set of equations are radiatively thick, 
coupled, or thin for a given flow region. 

The problem treated here uses as its starting point the problem 
treated by Bush [3], with the exception that radiation is included in 
the equations of motion, so that changes in the flow quantities in the 
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various flow regions due t'o radiation, as well as the need for addi- 
tional regions, must be considered. For this case the . Boltzmann num- 
ber (ratio of the convected energy to the radiated black body energy 
at stagnation conditions). Bo » (p U C 

CO OO 

the freestream Mach number, M = U^/Gyp^/p^) 2 , and the Reynolds number, 

R = (p^U go to infinity, and the density ratio, e - (y-1)/(y + 1) » 
the Bouguer number, Bu - o^a, and Q = Bur n+<,0 /eR (ratio of the molecular 
mean free path to the photon mean free path at shock layer conditions), 
go to zero . 

Since, consistent with the radiating hypersonic blunt body pro- 
blem, the freestream Bouguer number is very very small, it is required 
that a precursor field ahead of the body have a large scale thickness 

~ T J_0 r- 

of the order Bu - 10 body radii, Nelson and Goulard j,14j „ This 
extremely large region precludes any chance of applying a thin layer 
analysis to it. (A region such as this is physically consistent with 
the fact that an observer at a large distance from the body can "see" 
a body under these conditions due to the photons emitted from the 
shock layer of the body.) It should also be noted chat this "far- 
field precursor" is two-dimensional in nature and the equations are 
those of the field type. This is true because the addition of the 
equation of radiative transfer makes the entire set of equations of 
motion elliptic rather than hyperbolic, which is normally . true for 
the region ahead of the shock in hypersonic flow. In view of this 
extreme analytical difficulty, the far-field precursor will simply be 
assumed to be one-dimensional in which an exponent approximation is 


T )/aT 

Pl° ° 


is of order one. 
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valid. Thus the far- field precursor is found to have two regions 
which are designated as the exterior and the interior regions. 

As the flow approaches the body, more photons are absorbed by 
the gas and consequently the temperature, pressure and intensity in- 
crease. In the case under consideration, the gas temperature rises 
to the shock layer order of magnitude, although still smaller by a 
constant multiple less than one, Heaslet and Baldwin fllj . With this 
relatively high temperature, the gas is capable of absorbing and emit- 
ting in a thin region ahead of the shock structure. This region is 
referred to as the "near-field precursor" and is the first thin re- 
gion in the flow field, and thus the first that can be considered 
without using a one-dimensional approximation. This region has a 
scale thickness of order (i/Bur n )«l with the velocity components 
and density having only correction terms. Figure 2 gives the sche- 
matic for the flow region. 

The optically thin outer region of the shock structure, with 
thickness 0(1 /r), is expected to be the next region, but the in- 
creased temperature causes this layer to be thicker by f and con- 
sequently the corrections on the velocity components and density are 
different from the corrections for the conventional outer region. 

Since this region is now of the same scale thickness as the middle 
region, 0(r U) /R}«l, it is expected that they would have a different 
effective thickness. In the analysis it is shown that the middle 
region has an effective thickness equal to O(r cs) ,/R)log R/BuT n4tU . It 
should be noted that the outer region does not "see" freestream 
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infinity vhen looking upstream and thus the boundary conditions for 
the fluid mechanical shock are not found in the conventional vay. Now 
the precursor-induced conditions, due to the upstream heating of the 
gas, are the conditions to which the fluid mechanical shock must be 
matched . 

The middle region of the shock structure is the next region of 
flow after the outer region. The scalings and solutions are the same 
as those found by Bush [ 3 ] with the exception that a constant of 
integration appears in the temperature and pressure solutions. The 
equation of radiative transfer shows that this region, which has a 
scale thickness of Ofr^/R), is optically thin, Vincenti and Kruger 
M. (Since the inner region of the shock structure has a scale 
thickness equal to Ovef^/R), it will also be optically thin.) Since 
this middle region is optically thin, the intensity is scaled as a 
correction to the value it has previously obtained. 

The next region after the middle region is the inner region, as 
found by Bush [ 3 ] . The scalings, thickness and solutions of the re- 
gion are the same as for the nonradiatirig problem. The intensity is 
still scaled as a constant plus a correction. Since the inner region 
is thin in order of magnitude, it "sees’ 1 an infinity when looking down- 
stream, and therefore the fluid mechanical shock conditions, consistent 
with the flow conditions, are found as in the nonradiating case. In 
the strict mathematical sense, the inner region does not match to the 
radiation relaxation region, which is the next region in this analysis. 
In order to match these two regions, higher order terms must be in- 
cluded and the reader is referred to Bush M for the matching. 
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Between the shock structure and the shock layer a new region is 
present. This region is referred to as the radiation relaxation re- 
gion and is one which allows the intensity to equilibrate with its 
equilibrium conditions, Heaslet and Baldwin [n] . The fluid mechanics 
of this region are minimal. The scale thickness of this region is 
0(e/Buf 9 ), which is much larger than the middle region thickness, 
0(r W /R). This is a result of the fact that the photon mean free path 
is much larger than the molecular mean free path as shown by Bond, 
Watson and Welch [i] and Zel dovich and Raiser [22J . Also since 
BuT n »0(l), the following ordering 'sequence occurs: F^/R « e/BuF n «e. 

Since K = T^/eR -* 0 implies that the shock layer can only be inviscid, 
Bush [3j, i.e„, the possibility of the viscous shock layer, K « 0(1), 
is eliminated. The radiation relaxation region' "sees" an infinity 
when looking toward the body, and therefore a second set of shock con- 
ditions due to the radiation effect are found. Thus the near-field 
precursor region and the radiation relaxation can be thought of as the 
radiation shock, with its appropriate boundary conditions, and the 
outer region, middle region and inner region of the shock structure 
can be thought of as the fluid mechanical shock, with its appropriate 
boundary conditions, imbedded within the radiation shock. This double 
shock is referred to as a shock within a shock by Vincenti and Kruger 

[21] . 

The shock layer follows the radiation relaxation region, and for 
the flow conditions it is optically thick and in radiative equilibrium 
with the surroundings. The scalings and order of magnitude thickness 
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are the same as found by Bush [3] and the shock layer has two or three 
distinct layers for the same reasons as in the nonradiating problem. 
Table 1 summarizes the shock layer. The only difference is the in- 
tensity as an independent variable. 

TABLE 1 

k = r'VeR - o d = r w /e 5//2 R 

D = 0(l) inviscid shock layer 

+ viscous body layer 

D -* 0 inviscid shock layer 

+ inviscid body layer 
+ viscous boundary layer 


To summarize: the exterior region of the far-fiela precursor 

has thickness of 0(l/Bu) » 1 and the interior region of the far-field 
has scale thickness of 0(l/Bu) and an effective thickness of 
O(l/Bu)log » 0 (l/Bu) . The near-field precursor has a thickness 
of 0(l/Bur n ) « 1 and the outer region of the shock structure has a 
thickness of oO^/R) « 0 (l/BuT n ) . The middle region is oCr^/R) in 
scale thickness but has an effective thickness of o(r tU /R)log(R/Bur 114 ^) 
» 0 (r W /R). It is also true that 0(r“/R)log(R/Bur n ' h “) « o(l/BuT n ). 
For the inner region the thickness is oCer^/R) « ©(r^/H) since e -• 0. 
In the relaxation region the thickness is 0(l/Buf n ) » OCr^/R). 


JPL Technical Memorandum 33-687 


T 



Chapter II 


THE EQUATIONS OF MOTION 

Consider the hypersonic flow of a viscous, radiating gas past an 
axially symmetric blunt body. An orthogonal curvilinear coordinate 
system, with the body surface as reference, is used. Let - ax and 
= ay, respectively,, represent the distance measured along the body 
from the stagnation point and the distance measured normal to the body 
at the corresponding location. The length, a, is the nose radius 
of the body. Figure 3« 

With the above non-dimensionalization, the curvature of the body 

surface and the body meridian radius become, respectively, k = «/a 

and = aB. The velocity components in the x^ and y directions, 

respectively, are u = U u and v = U v and the pressure, density, 

temperature, and intensity, respectively, are p = p p, p = p p, 

X °° X co 

4 4 

T ± = TJT, and L, = crT Q I. The quantities U^, p^, p^, T^, and oT , 
respectively, are the free stream velocity, pressure, density, temper- 
ature, and the black body radiation present in the stagnation region. 

A perfect gas is assumed (p = pT) 'having constant specific heats 

(C j C , and y = C /C = const.). The Prandtl number is assumed 
1 P 1 P 1 V 1 

to be of order one, P = 0(1), with the viscosity coefficient being 
proportional to the power, cu, of the absolute temperature, 

(w>^ * = t where j £ cu £ 3/2. The radiative absorption coef- 

ficient is assumed to be independent of spectral frequency (the gray 
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gas assumption) and proportional to the first power of the density and 
the power, n, of the temperature (o' = a a = a pT n ). For molecular 

1 CD CD 

radiation of a monotonic gas, n * 6, Thomas [is] , Goulard [V], and 
Penner and 01 fe [16], 

The Navier-Stokes equations are assumed to be valid. The radi- 
ative pressure terms in the momentum and energy equations have been 
neglected, since they are negligible in all regions of flow considered, 
Goulard ^8^. The addition of the radiative term to the energy equation 
is therefore the only change that is necessary to the continuity, mo- 
mentum and energy equations that are valid for nonradiating flows. An 
additional variable, the intensity, requires the inclusion of a radi- 
ative equation of transfer. This is accomplished by using the Milne- 
Eddington differential approximation: 
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The components of the radiative heat flux vector are eliminated by 
differentiation in favor of the integrated intensity, thus yielding: 



Traugott, [19] and [20J , performed a shock structure analysis by 
eliminating the intensity in favor of the heat flux vector, which 
had been assumed to be one-dimensional. The equations of motion 
under the above conditions are; 


Mflll + I 5 (p u l + £21 + 

By h dx h 


linqj + vcos'p) 


- 0 , 


du , kuv' 

+ V— C — + -7— 

dy h 


dv usincp 1- vcoscp* 
dy r 


2s incp 


( u 3v , dv 1-e 1 dp 

h dx - v dy " h / 1+e ^ dy 
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usirup + vcoscp 
r 



oBu 
B 


(i - 4r"V*) 


> 


1 

h 



JL — 

r a byj 

\ 2 
coscp ] = 3 <*Bu. 


( 2 . 5 ) 

(i - 4r“V) , 


P = pT , 

( 2 . 6 ) 

_U) 

H = T , 

( 2 . 7 ) 

a = pT n j, 

( 2 . 8 ) 


JPL Technical Memorandum 33-687 


11 



where h = l + *(Y + y),r = B + (Y + y) coscp, e= y-l/y+1 - o (the 
Newtonian approximation), m 2 - evJ/y p. » 1, R = P^U^a/jjf » 1, 
the Boltzmann number, Bo = p U C T^/aT ^ = O(l), T = T /T >> 1, 

r oo oo p O' o '/•» o » * 

and the characteristic optical length, Bu" 1 » 1. 

Along with the above equations of motion one additional equation, 
in the form of a boundary condition, must be satisfied in conjunction 
with equation (2.5)- This equation is: 


SI 

dy 


lim 


= /3 Bua (I. - Ur'Sh,* 1 ) 

' lim B 


(2-9) 


where L. . and 

lim dy 


lim 


are taken in the limit as they approach the 


downstream boundary of the regiun, and T is the temperature at the 

J3 

upstream limit of the region immediately downstream. Equation (£-9) 
is implicitly satisfied at each boundary. 


12 
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Chapter III 

THE EXTERIOR REGION OF THE FAR-FIELD PRECURSOR 


With a body in hypersonic flow, traveling sufficiently fast so as 
to induce electro-magnetic radiation, the flow quantities start devia- 
ting from their freestream values prior to their encounter with the 
"shock wave." The precursor effect is due to the generation of 
electro-magnetic radiation in the shock layer, which is propagated 
upstream at its characteristic speed of light. 

The first region of the far-field precursor is the exterior 
region. The scale thickness of this region has an order of magnitude 
thickness, (Bu)~\ equal to 0{l0‘ 1 ' 0 ) body radii. Although this region 
is two-dimensional, it will be assumed that it is one-dimensional in 
order that analytic solutions may be obtained. This assumption is 
valid since the changes in the flow quantities are small and therefore 
would not significantly affect regions that follow downstream. This 
region has the pressure and temperature varying as corrections to 
their freestream values. The corrections to the velocity components 
and density are very small but are included for completeness. The 
intensity is scaled from its freestream value. The coordinates and 
flow quantities consistent with the above description are: 


y = B^ V 


-3 


er 

u = coscp + " u + — - 

w e 

r “ 3 

v - -sincp + -5- v + — 
n e 


(3.1) 
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3 


(3.2) 


p - i + r p + — • 
e 


0 ' 1 * 3" °e + — ' 


T»i+r J T + — , 
e * 


i = Ur i + — 
e 


Substitution of (3.1) and (3.2) into the equations of motion, 
(2.1)- (2,8), yields, to first order, the following: 


P e - T e ’ V e " p e sint P = 0 


dv 3u 

" Sinfp ^e + ^ = ° 5 - sinrp ^T 


y r-2- = c 
e by 


- sla * SjJ - B? (h - v 


- 3 (I - 1) 


The solutions to equations (3.3) are: 


p = T = e“ /3 y e 0 

e /3 Bosirxp 3 "e 


— T T e » 
sin cd 


sincp e e sin rp e 


I = i + J 6 


~/3 v 
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where it should be noted, due to the one- dimensional assumption, that 
J is not a function of x. 
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Chapter IV 

THE INTERIOR REGION OF THE FAR-FIELD PRECURSOR 

The next region is the interior region of the far -field precursor. 
In this interior region, the intensity is much larger than its free- 
stream order of magnitude, but still small compared to the stagnation 
order of magnitude it attains downstream. The pressure and tempera- 
ture are still of their freestream order of magnitude and are scaled 
correspondingly. The velocity components and density are again scaled 
as the fluid mechanical aspects of this region are still small. The 
scale thickness of this interior region is the same as that of the 
mctGYiov regionj 0(Bu )> tut has an effective 'thickness order 

of magnitude which is much larger. Hence the coordinates and flow 
quantities have the following representations ; 

y = Bu y J J 

U = COS<p + ^ U J + , 

v = -sincp + — + — J 

P + Pj + > (b.2) 


1 6 
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T = T, + — , 


i = 4r i. + — . 
j 


Substitutions of (4.l) and (4.2) into the equations of motion 
(2.l)-(2.8), yields, to the leading order, the following: 

p = T , v - p sintp = 0 , 

J J J J 


dv . dp. du 

- sin-p r — ^ + r — ~ = 0 , -sinm -r — ^ 

ay j dy j Sy j 


9p 


- Y 


J ay, 

J 


= o , 


ST. 


-sincp 


n 


ay, ' Bo "j s 

J 


d 

ay 




n 


= 3 T I. 

3 A A 


(4-3) 


The solutions are most conveniently found by transforming the inde- 
pendent variables from y. to T.. Hence the solutions, in terms of I., 

0 J J 

ares 

= t j * 




. 2 5 
sm cp 


(4.1.) 


V. 

J 


sincp j 


T 4 , 
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Y’ 

U. H 2 - T. , 

J s incp j ■ 


, / 3 -B o ^i ny . ( , 1} _ 


0 X • j 

T. 




A. 

J 


4j. 
i 


where, for convenience, A. = 1 + w v “ ; — . 

J / 3 Bo sincp 

The solutions in the limit as y -» <» are; 

J 

4j. ^ 

p = t = 1 + 3 y i 

j /3 Bo sin® J * 


p -i 2 » 


T. 


sin (p 


V = — i — T 
j sin? j * 


Y\ 


u . = 


s in® j 


V > 


(*. ; 


^ ^ = V3 ^ . 


and in the limit as y . -* 0 : 

0 


Pj ’ Tj " V n y.) 1 '” ’ 

J 


V = - - T 

j sin® j 5 


(4.6) 
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*<n 


1 




sm cp 


- 


U. - - 


Y\ 

ii , 

sirup j 


/*3 Bo sincp 


(/3 n y^) 


1/n 


The matching of the interior region of the far-field precursor 

to the exterior region of the far-field precursor is accomplished 

using the normal coordinate as the matching variable, i.e., 

= (K.Z. - y)A n where Bu X "* Through the matching, it is 
J J -L J- 

found that J, - J (not functions of x), K. - l//3 and the effective 
J ® J 

thickness, Z^ = (l/Bu) log 
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Chapter V 

THE NEAR-FIELD PRECURSOR REGION 


In the near-field precursor region, the intensity and temperature 
rise to their shock layer orders of magnitude, vhereas the pressure 
has risen to the order of magnitude it has in the middle region of 
the shock structure. The thickness of this region, 0(l/BuT n )‘, is 
much less than one, but is much larger than the order of magnitude 
thickness of the shock structure. The velocity components and 
density vary as corrections to their freestream values. Hence, the 
coordinates and flow quantities have the following representations : 


X - x. r , y - <5Y r (x r ) + — y T 


BuT 


u 


COSCp + € u r + f 


(5-1) 


v = -sincp + ev^ + 


p * rp r + — , (5.2) 

p = 1 + ep r + — - , 

T * TT r + — , 

1 = hi + — . 

T 
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Substitution of (5*1) and (5*2) into the equations of motion, 


(2.l)-(2.8), yields, to the leading order, the following} 

v r " P r sin<v = °' P r = T r * 



(5-3) 


.The energy equation retains the convection and radiation flux diver- 
gence terms with the pressure gradient, heat conduction, and viscous 
dissipation terms being of higher order. 

The solutions to equations (5*3)> in particular the energy and 
transfer equations, are difficult to find due to their nonlinearity. 
Although they are nonlinear, they can be solved since they can be 
transformed into a modified Bessel equation. By letting 
d/dT| = (l/T r n ) d/3y r and combining the two equations one obtains: 

- sint p s# = 35s If ■ (5J * : 

Integrating once and transforming from 7] to T r as the independent 
variable, by use of the energy equation, yields: 
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(5-5) 


k SI 

K T = (I - T 4 ) —■ 
r v r r ' 3T 

r 


2 2 2 

where K = 3Bo sin y/l6 . Letting Z = the equation becomes 


dz = - (I - Z 2 ) 
dl K {I r Z } 

r 


(5.6) 


which is the Riccati equation. The Riccati equation can be trans-. 
formed into the Airy equation (a modified Bessel equation) by letting 


Z = (K/2) u’/ u j where u' = du/dl^ . Thus equation (5*6) becomes: 


u" - — I u « 0 . 
K 2 r 


(5.7) 


The solution of the Airy equation is in terms uf modified Bessel func- 
tions. Obtaining the solution of (5*7) and transforming bach into 
temperature yields, after imposing the condition that -* 0 as 

X - 0: 


T r = 


i ihd 


J*_ t 3/2 
3K r „ 

. jjtj m 

-1/3 3K r 


(5.8) 


where the numerically subscripted I is a modified Bessel function 
and the letter subscripted X is the intensity. 

Therefore the general solutions for the near precursor region, 


in terms of 1^, are: 
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T = 
r 


i 2/3 


JL I 3/2 
■ 3 K r J 


I riL I 3/2] 
I -l/3 3K r 


p = T , 
*r r ’ 


v = — F- T , 
r s incp r 


— - — T 
. 2 r 
sin cp 


2 Y' 

-r-^ T , 
sinrp r 


r 

>inp / dr 

1 1 

J T (I “ ^ 


The solutions in the limit as -* » are: 


p = T = 

* r r 


r r 


C/3 n y r ) 


T7S + ’~ ' 


v * — F— T , 
r sincp r 


— i — T 

. 2 r 
sin 9 


T , 


2 Y 1 

— r— T , 

sincp r 


I _ /~3 Bo sir 


C/3 n y r ) 
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0 are: 


and in the limit as v 

r 

« t * e (x ) 

*r r r v r 


ho* 


r Bo sincp v "r 7 1 r 


(J “ ©v. ) y« + — 


V = — — rn 

r sincp r * 


= 


T 


. 2 ~r * 
sm <p 


( 5 - 11 ) 


u = - 


a r 


T 


r sincp r * 


J r - W + b r <* r > » r “ y r + — 


Applying the boundary conditions at - 0. the following re- 
lations for J and b in terras of 0 were found: 
r r r 


<v- 


, * h/3 

it- J 3/2 

J r. J 

r T 

L J 3/2 
3K r 


( 5 . 12 ) 


and 


b s . ig o sincp . 
r X ®r ‘ 


The matching of the near-field precursor region to the interior 
region of the far-field precursor is done with respect to the nor mal 

st *- n 

coordinate, i.e., y g = (- eY m + y)/x g where BuXg -* 0 and 
BuT n X g -* ». The matchings are found to he identical, and thus 
no new constants are determined. 


2h 
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Numerical results of this region are also presented in Figure 4, 
in the form of a plot, along with. the numerical results for the radi- 
ation relaxation region, to be discussed later. In this figure the 
fluid mechanical shock wave is represented by a discontinuity. These 
results can be compared to the results of Heaslet and Baldwin [11], 
for e 'f 0, given in Figure 1. 
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Chapter VI 

THE OUTER REGION OF THE SHOCK STRUCTURE 


The next region downstream from the near-field precursor region 
is the outer region of the shock structure. The scale thickness of 
this region is 0 (T^/r). This scale thickness is larger, by a factor 
r* 0 , than the conventional outer region, since the temperature has 
risen to its shock layer order of magnitude in the near-field pre- 
cursor due to the effectiveness of radiative transfer. In this re- 
gion, all variables are scaled as constants plus small corrections. 
The fact that the pressure and temperature are varying as corrections 
to constant values, as opposed to variation as free variables, is the 
essential difference between this outer region and the conventional 
outer region. The coordinates and flow quantities, therefore, have 
the following representations: 

X = x o , jr - e Y 0 (x 0 ) + §- y 0 5 (6.1) 


u = coscp + 


n l/p 
Q ' u 






v = -sincp + f 

J> - r |n o (x o ) .ap o ) + — , 

P - 1 + 0 3/4P p Q + — , (6.2) 
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t = r 

i 


(e o (x o ) ♦ n T o j 

- 11 ( J o<*o> + n Z o) 


+ 


where fl = BuT n+U) /R 

Substitution of (6.1) and (6.2) into the equations of motion, 
(2.l)-(2.8), yields to the leading order the following: 


n o = 9 o ' 


p = T , 

•o o * 


- p simp = 0 , 


'•sincp 


dv 

W, 


o ,1-3 A* dP o _ k „U) ’ o 


d 2 v_ 


+ 2e Q 


5y o = ^ ®o 2 


du 




. ™o -a^S, 

-sincp 2e Q ' Y ' r — = 0 ~ 

5y o ° ay o °ay 2 

o 


r l Y;cn iAp , 

y o 

. 3T o ®o A o. + 4 Q n ,, .4, 

-° lr “P »T - F — 5 + 5S 0 o (J o ‘ ®o ) 

° Sy„ 


< ay o 


The solutions to (6.3) are: 


(6.3) 
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p = T ® S (x ) A exp 

*0 o o v o' o * 


4 e n (j - 9 **) 

A ' A A ' 


P since 
(U y o 
. 0 o - 


. . , 1 * e n4 ”(J - e ■*) 

o' o o' . o o o' 
Bo smep o 


Bo P Bln tp 


v ° ' A ° exp [■ T? y °] 

L o J 


+ e n l- 3 AP f £!° 


lfP/3 


a ~~ f j r 

> o _ I Psirup 

L* 9 “ 0 


8 e n (j -eS 

u u o 

~ ; — . 2 — y 0 

Bo sin cp 


8 fi n " H “(.T - e *) t lf 
- °3 ° <!-?>> 


Bo sin~cp 


P„ - 


o sincp 


”o ' v o (x ^ A o k/3 eJ ®|- y o 

9 o 


{• 


(6A) 


- e n 1 ^ X' A 


o o 


exp 


[ . iglag v 

^ e o °j 


- € 


n 1_1 / p | 


2 S o A o * P /3 

Bincp(l-P ) eX ^ > 


Psincp 

e 10 y ° 

o 


6 e n (j - e 6 e n+u) (J - e **) 

o o o' + ^0^0 o ' 

2 3 

Bo sin tp Bo Bln (p 


(1 + 


|)J , 


I = a + "b 0 n y 
o o o o •'o 
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The matching of the outer region of the shock structure to the 
near-field precursor is performed using the normal coordinate as the 
matching variable, i.e., y^ = (-eY 0 + y)A^ vhere and 

(r'VR)!^ ”* “• During the course of the matchings it is found that 

e o (x o ) ■ e r (x r ) ’ W = J r (x r ) 


and 


fc.OO . b.(o - - ^ 0, 


( 6 . 5 ) 


It should be kept in mind that J and 0 are not independent con- 

r r 

Btants. It was also found, through the tangential velocity matching, 
that Y^ = 4/3 ^ • This is interesting, since this result implies 
that the near-field precursor region decreases in size as one pro- 
ceeds away from the stagnation streamline. Stated more simply, the 
near-field precursor region is a radiation "cap" that sits ahead of 
the shock wave and becomes smaller in relative size as one moves 
along the shock away from the stagnation streamline and then finally 
disappears. Of course, the present theory breaks down at some point 
prior to the disappearance of the near-field precursor, since it is 
assumed that this region is to be of significant size. 
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Chapter VII 

THE MIDDLE REGION OF THE SHOCK STRUCTURE 


The next region downstream from the outer region of the shock 
structure is the middle region, or dissipation region, of the shock 
structure. The scale thickness of this middle region is O^^/r), 
which is the same as for the outer region, but has an effective 
thickness which is much larger. The pressure and temperature are 
now free to adjust upwards, due to the fluid mechanical aspects of 
the shock wave. The velocity components and density are still of 


.their freestream order uf magnitude. The intensity is still scaled 


as a constant plus a correction, since the region is optically thin, 
which excludes any change of order one or greater. Thus, the coordi- 
nates and flow quantities, which (except for the intensity) are the 
same as for the nonradiating problem, have the following represen- 
tations: 


r WJ 

y = + tt 


R y m 


u « u + — ■ , 

m 3 

V = v m + > 

P = r P m t , 



(7.2) 
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t = nr + -— , 

m 


1 = 4 


[J (x 
l m' j 


) + m: 

m' R 


n-Hu 


m 


) 




Substitution of (7*1) and (7*2) into the equations of motion, 
( 2 . 1 ) -( 2 . 8 ), yields to the leading order the following: 


p = p T , 
m mm 


Pm v m = - Slnt P > 


v sin-p + -=■ T 5 — = -sin t , 

m r 3 m ay r * 

m 


(7.3) 


du 

u^sincp + g — = coscp sincp , 

y m 


ST 

m , 

sincp v — + 


1 - ( T w &Tjn V T m - (* Vm \ ! ( 6 rn \ 

Pay m\ md ^m/ m 3 \ S W \ dy m/ 


* 0 , 


d y, 


m 


x K m m r 


Equations (7*3) are most easily solved by using a variable transfor- 
mation. A modified Crocco transformation is used in which the inde- 
pendent variables (x^, y^) are changed to (x^, V^) where 


V = v + since - 
m m 


(7*4) 


Equations (7*3) in terms of the new variables (x m , V^) are: 
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> 


P = P T 
ra m 


p v = -since , 
mm r ’ 


SU m k U m 

m 4- m tt «_ 

TT? T7- = 0 , u = u ~ cosrp s 

oV w 3 5 mm * 

m m 


S 2 T P - 1 ) 5 T 

m '3 n 


fV „ 4 / dU mV 

av p 3 + l av j 

m \ m / 


- 0 , 


*» 4 T 


-ZJ2 + i J2 =0 
av 3 v u 5 

m m 


a / V 51 
o / m m 

5V_ l m n+tn 5V 

ni \ u i m 

\ r m m 


) =° 


The solutions to ( 7 * 5 ) are: 


sincp 

v 3 
m 


u ~ coscp + w (x ) sincp + v 


T - OJxJ + S (x ) FsirKp + v ] 4p / 3 
rn m m m N nr ^ ~ mj 


p r . 1 2 

3^2P L SinCp + V mJ 


m r . 1 8/3 

sincp +v ' 




y = L » T - 7 
m m 3 sincp 


• sincp + A 
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n+u) 


4/b (x ) 

( ) m m 
m ~ nr nr ’ 3simp 


m 

f A_ 

I sin .9 + 


-sincp + A m 

where now 0 m appears in the solution to the temperature . It is 
needed since the gas is preheated due to the precursor effect. 

The matching of this middle region to the outer region of the 
shock structure is most conveniently done with respect to the normal 
coordinate, i.e., y^ « (-eY^-K^Z^+y) X ^ where (r/t <W )\^ "♦ ® and 
(I0'*)K Z » (R/r W )Xt,- It is also required that (eR/r^^U -o so 
that the normal and tangential velocity components will match for 

To perform the matchings, the solutions of equations (7.6) in 
terms of (x ,y ) as y -* (R/r tW )K Z -* oo are found to be: 


since - A (x ) exp 
T mm 


y 

u) 


u = eosep + w A exp - y s 

m Y mm ato m 


T - 0 + S A exp - Psi ”^ y , 

m m mm qud m 

m 


p = p T , 
-^m m ' 


v - -sirup + A exp - y m , 

m t m * i - (u v m 


I = a + b 6 
m m m m m 
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Performing the matchings it is found that 


j(x)=j(x)=j(x). 
m m ' o' o' r r i 


©JO = = Q ( x ) 


m' m 


o' o r r 


S n (x m ) “ s o (x o> > V x m > - A o (X o J > 
*> CO - b Jx) = . a ° £ l sa> e. . 


in nr 


o' o' 


( 7 . 8 ) 


It is also found that the restriction for matching on* the Prandtl 
number of 5 < P < ~ is required as it was for the nonradiating pro- 
blem. The effective thickness is found to be Z ffi = (r^/Rjlog R/Bur n+U) 
with = O^V'Psincp. 
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Chapter VXXI 

THE INNER REGION OF THE SHOCK STRUCTURE 

The next region of the shock structure is the inner region which 
is thinner than the middle region by a factor e. The tangential velo- 
city component and temperature have already reached their stagnation 
orders of magnitude and are therefore scaled as corrections to their 
middle region values. The density finally rises to its stagnation 
order of magnitude, which must be accompanied by a corresponding rise 
in pressure and a decrease in the normal velocity component. The in- 
tensity is again scaled as constant plus a small correction, since 
this region is also optically thin. The coordinates and flow quan- 
tities in this inner region have the following forms; 

er m 

x - x t , y = eY i (x i ) + — y i (8.l) 

u = W. (x. ) + eu. + , 

v = ev. + , 

t ’ 

P = (r/e)p. + — , 

P = (VOpj. + S (8.2) 
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T = r^fx. ) + eT.j 
I - U (j. ( X .) + 0 


+ --- 


+ -«*» 


The equations of motion (2.l)-(2.8) consistent with the above 


expansions are: 


p i = p i 9 i * 


P i V 1 = -sincp , 


o • 4 A u) dV i . 2 

-2smcp ~ - - 0 i = sin cp , 


e U) 1 M / _ _ x 

■■■■ — ZZ q rvn f ^ X‘J ! 


@ w BT, r n 

|r- = sin «P (|sin 2 tp - G i + 0 m ) + ? (coscp - W ± ) S J 


1 SI i 

Pi% n Sy * 


= 0 . 


vhere - v 


p i “ 


Y ! W. . The solutions to these equations are : 


sincp 

V, ' 


p i “ P A 


= w j( x j) " sincp (coscp “ W^) 


{Y i (x i ) - y ± \ 
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+ Kcos:p - W.) 2 J- Yi 


e 


/ 20 . \ 

. = a . (x . ) + rr— r (lOg I V. + —7 — - j ) , 

i i i' 3sincp V 1 sin( P / 

r . 26 . 

p simp , + 

p 20. '“ i sincp ; 


sin (0 (V. + — r~) I = 
v i srncp ' 


3ii2-2 \ x 

o 0 1-nu * 


i - y il] 


The matching of the inner region to the middle region is most 
easily accomplished by using the normal velocity as the matching 

— *• — — t» 

variable, i.e., = v/£^ where 3^-0 and -T./e - ®. In performing 

the matchings, it is found that: 


J. (x.) = J (x ) = J (x ) , 

i' i' rrr nr r v r' 5 


and 


b. (x. ) 

i' i' 


a (x ) 
m nr 


w (x ) 
m m' 


S (x ) 
m m 


- ejxj 


m' m' 


-e 


r<*r> f 

.sinrp J 


kb 

3sinrp 


-sincp + A 


r • T - 


P T 

m m 


n+u) 


dv 


sintp + v m > 
m 


(cos:p - W.) 

(sinf) 4 ^ * 

8j - 8„ * 3^p s 1 ” 2 * + 2x£fi Ccos ? ~ V i> 

(sincp) 4 ^ 


( 8 . 
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Chapter IX 

THE RADIATION RELAXATION REGION 

The region that follows the inner region of the shock structure 
is the radiation relaxation region. This region was not present in 
the nonradiating problem, since there was no need for a region in 
which radiative equilibration would occur. The rise of radiative 
intensity in this region is accompanied by a decrease of temperature, 
Heaslet and Baldwin [ll] . This decrease in temperature, or cooling 
of the gas, is referred to as radiative cooling. This transfer of 
energy, in the form of temperature, out of this region is responsible 
for the precursor temperature rise before the shock wave. The flow- 
quantities in this region have all reached their shock layer order of 
magnitude, and are therefore scaled as such. The scale thickness of 
this relaxation region is of the order of the photon mean free path 
at shock layer conditions and is thus equal to o(e/Bur n ). As this 
region is mainly a radiation- dominated one, it would be expected that 
all three terms of the transfer equation are presentj i.e„, gradient, 
absorption and emission terms. This is, in fact, the case, as will 
be seen in the equations . 

The coordinates and flow quantities which describe this radiation 
relaxation region are; 
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* ■ x a 


y ’ £Y a (x a ) * — H 

Bar 


(9-1) 


„ 2 

u - W + e u + 
a a 

v = ev + , 

a ’ 



(9.2) 


i - r t + — , 

a 

I = 4l + — . 

a 


Substitution of equations (9«l) and (9.2) into the equations of 
motion, (2,l)-(2.8), yields: 


<0, V J - K 4; (p A> = 0 > 


du dp 

- ElKp 5^r - »a sr * 0 > 


dv dp 

a 

-sincp t- — + 2 r — - 0 

5y a & y a 


n a - Pa T a • 


(9.3) 


st 

-sincp 

a 


a ipT fl (l -tS 

y. Bo p a a v a A a ' * 


a y c 


W. B *v 


3 p a T a" * .O 
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The solution to these equations was carried out numerically and a 

representative solution is shown in Figure 4, as mentioned before. 

In order that upstream and downstream matchings might be obtained, 
the solutions in these limits must be found. The solutions in the up- 
stream limit, as v -* y are: 

7 a a 7 


I 

a 


= J (x ) + b (x ) 
a a a a 



e 


n-l 

a 


(y. 


V 


9 


, 0 . 4Sina e 

a a Bo i 


n-l 


( J a‘ 


O (jr a 


V 


p a = I (c x - v a sin<p) , 

_tk 

T 7 (9*4-) 

cL 


u 



2Y T p ) 
a a 7 


9 


V 

a 



+ Y* W 
a a 


and the solutions in the dcvnstream limits as y -* are: 


X a * (x a } ‘ B(x a> ex T ^ <2 9 B “" 1)y . ' 


T = 9 + T — ^ 

a B /3 Bo sincp 


exp 


/3 




5 


P a = i ( c ! - v a sinc p) 


(9«5) 





> 
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1 

sintp 


i 


u a = “ 


(= 2 + a; p a ) 


V 

a 


smcp 

P« 


+ Y* W 
a a 


The matching of the radiation relaxation region to the inner re- 
gion of the shock structure is accomplished using the normal coordi- 
nate as the matching variable, i.e., y^ - (eY & - y) X ^ where 
(EAr“)X 6 -* co and (Buf n /e)Xg "* 0. The matching establishes the up- 
stream boundary conditions for the radiation relaxation region, and 
therefore this matching finds the ’'fluid-mechanical*' shock conditions, 
(it should be noted, as mentioned before, that strict matching for 
these two regions was not accomplished since higher order terms must 
be included. The reader is again referred to Bush for this match- 
ing.) Through performing these matchings the following relationships 
are found: 


0 (x ) = 0. (x. ) = 
a v a 7 i v i / 2 m 9 


\(x a ) = W i (x i ) = coscp , 

n a (x a ) = s ( x 0 ) = J ,( x i) = j J x J > 


( 9 . 6 ) 


a' a 


i' x J 


r' r' 


and 


W ■ W - - iS r fla 8, • 


An interesting result is that the pressure is constant, to leading 
order, in this region. This fact was first .realized by Goulard |loJ , 
and later discussed by Penner and Olfe , It should also be noted 
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that the tangential velocity and pressure have attained their Rankine- 

Hugoniot conditions at the start of this region, whereas the tempera- 

ture is above the Rankine-Hugoniot value of sin cp/2 by 0 m « Since 0^ 

is also equal to the amount of temperature rise before the shock, due 

to precursor radiation, it can be seen that the temperature jump 

across the fluid mechanical shock only is Rankine-Hugoniot to the 

leading order. In other words, the temperature jump across the shock 

is sin cp/2, i.e., from 0^ to 0^ + sin cp/2. In the present analysis 

the equations of motion in the radiation shock structure and fluid 

mechanical shock structure are one-dimensional, and consequently the 

conservation of mass, momentum, and energy must be maintained in the 

normal direction. This fact means that the jump conditions across 

the radiation shock (with the fluid mechanical shock imbedded within} 

must also be Rankine-Hugoniot, i.e., a temperature jump from 0 to 

sin cp/2. Therefore the downstream value for the temperature, 0 , is 

B 

given by: 


0 B = 0. 


9 

°m 2 


( 9 . 7 ) 


In this analysis the radiation and fluid mechanical shock struc- 
tures were one-dimensional, as mentioned above, as a consequence of 
the assumption that I/BuT 11 « 1* If, on the other hand, the assump- 
tion of l/buF n = 0(1) were made, the radiation structure would becane 
two-dimensional,, Two-dimensionality in the radiation structure allows 
transfer of energy in the longitudinal direction® Therefore, conserva- 
tion of energy in the normal direction is therefore no longer a 


U2 
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requirement, and this manifests itself in that the Rankine-Hugoniot 
conditions across the radiation structure are no longer valid. The 
jump conditions across the radiation structure are less severe vith 
energy loss in the longitudinal direction. It should be noted from 
physical considerations that the fluid mechanical shock structure, 
imbedded within the radiation shock structure, is still one- 
dimensional under the assumption of l/Bur n = 0(l) , Therefore, 
Rankine-Hugoniot jump conditions across the imbedded fluid mechanical 
shock structure apply locally in both the one-dimensional and two- 
dimensional cases, although the initial and final conditions are 
different. In the two-dimensional case, the radiation relaxation re- 
gion thickness is 0(e/BuT n ) = 0(e), which is equivalent to the 
thickness of the shock layer, and therefore the radiation relaxation 
region extends all the way to the body and has terms common to both 
the radiation relaxation region and shock l&yer in the one-dimensional 
case. 
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Chapter X 
THE SHOCK LAYER 


The region between the radiation relaxation region and the body 
is the shock layer ana it is inviscid, since K = r^/sR - 0. In the 
case under consideration, the shock layer is in radiative equilibrium 
and the orders of magnitude as determined by Bush pij, Chester , 
and Freeman £f] are appropriate. Thus the coordinates and flow quan- 
tities. have the following forms; 


x = x T 


( 10 . 1 ) 


u = a L + > v = ev L + > 

p - r pl ♦ — > " " I p l + — 

t = r *r + — , i = 4i^ + — , 


( 10 . 2 ) 


The leading terms in the equations of motion (2.1)- (2.8), con- 
sistent with the above expansions, are: 


Pl = P L T L 




sir »l * sir (B “h V ' 0 • 

b li 

2 S^- “'’Lb . 2 ’ 0 ’ (10 - 3 ) 
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yields tne following solutions for the inviseid shock layer ; 
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+ _ 1 /, + 2 2v4 

h is h - h a ) 


The matching of the shock layer to the radiation relaxation region 

-o — t 

was carried using the matching variable y^. = (eY^ - y)/\y where 
(BuT /e)X^ - ® and x^/e - 0. These matchings give the "radiative" 
shock conditions. Thus the shock layer boundary conditions were 
found to be: 


“t ( V V = ”a (x a ) “ COS, f ’ 


*t <*l« y l> ’ nt-J - #2 . 


<*L> r L> ‘ 1 


. 2 


(10.5) 


T r x y ) = sin 2 

T L ' V 1/ 2 J 


v L (x , Y ) = Y t * coscp - sin^j 


. 8 

T_ (x Y ) = 

T, V L> L' ^5^ 


h6 
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Chapter XI 
THE BODY LAYER 


The coordinates and expansion appropriate for the body layer are't 


x = x 

c 



( 11 . 1 ) 



c 







i ” r i + — 

c 


( 11 . 2 ) 


I = 4l 


The equations of motion consistent with the above expansions are: 


p = p T 


* 



< B P c u c> 



(B p v ) = 0 
K c c 1 


9 



(11.3) 
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3u 3 a 


F VW -v ' V — 

c c c 

c c 


dp ^ 

/ c C\ ^ c ^ o /„ u> C\ 

P ( % — + v - — -) + 2 ^ — = D r — (T t — ) , 

~n ' n Av 4v ' Air ' r> A.\r * * 


du 


dy c c 3y c 


3T 3T _ . 3T 

p ( u _£. + v _£) = H _L_ ( T « _£) 
c ^ c 3x c 5y ' P 5y ' c Sy ' ? 

c c c c 


X = T , 
c c i 


where 


ja 


D = 




The boundary conditions, obtained from the matching of the body 
layer to the inviscid shock layer, in Crocco variables (s^ = x^. 
t = u /coscp, t = |T 7b[ {du /by | ), 


are : 


P c (Sp^c - "0 * P„(0 


C' c' 


T c (“c^c “ 1 i > 


T (s ,t -* ») , 2 ^ ^ P (s ) 
c c c 


c ' c' 


(11,4) 


I (s ,t -* co) = 
c v c’ c ' 


TE 


The solutions to equations (11.3) for the inviscid case, D -* 0, 
consistent with boundary conditions (11.4) are: 

P c (» c .t c ) = P c (s c ) . 


T (s ,t ) = f , 
c v c* c / 2 ? 
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J 


T 

C 


(s ,t ) 


= 2 1_iJJ P (s ) 
c c 


f B + \ _ COSCP 

' c'V 2B(s c )p^(s c ) 


t 



x c (s c’ t c ) " h ’ 

t (s ) = [-2 ln{2 P c {s 0 )} . 

C w c 7 -k— - 

* n. osfrt 


(U.5) 


The boundary conditions for the viscous body layer are satisfied 
in the same way as done by Bush £ 3 ] . However, since the viscous body- 
layer solutions are numerical, they will not be presented here. 
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Chapter XII 

THE VISCOUS BOUNDARY LAYER 


The coordinates and flow quantities for the viscous boundary- 


layer are: 


X=X BL > 



( 12 . 1 ) 



( 12 . 2 ) 


P 


1 

e P BL 




> 


T=rT BL *-~ , 
I = %L + -~ • 


The equations of motion resulting from the above expansions are: 


P BL P BL T BL 


s 


a 

SX BL 





(B P BL V BL ) *° 


Sp. 


BL 


3y. 


BL 




( 12 . 3 ) 
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du. 


BL 


P BL Wl dx^ ' 'BL dy. 


+ 


* V, ^£BL _ _d_ , to dU BL v 

) + ^ Av “ Air ^PT Air I > 


BL 


SX BL ay BL Bl dy 


ST. 


ST, 


3L 9^ BL Sy BI / P 8y BL ' BL 8y. 


BL 


BL 


rn 

/ ‘'BL ^ ^ J ‘BL\ _ 1 _d_ a) X 

PfiL ^BI ^' w ' ^tr l P Sv V^PT Air / * 


^L 


= T. 


BL 


The boundary conditions at the outer edge of the boundary layer, 
obtained from matching with the inviscid body layer, are: 


BL 


(x BL’ y ] 


BL 


°°) ~ i > 


^BL ^ X BL , ^BL " ®) ~ ^L^L^ ~ P c^ X c^ * 
“BL (x BL< y BL " -) = [- 2 ln . 

•^L ^BL^BL "* = T6 * 


(12.4) 
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APPENDIX 


Numerical Solution of the Radiation Shock Structure for 
Four Values of the Boltzmann Number 


The energy and transfer equations for the near-field precursor 
and the radiation relaxation region are solved numerically for a range 
of Boltzmann numbers . The equations to be solved in the near-field 
precursor, see equations (5*3)> are: 


d T 


-sincp ~ T 0 (I - x ) , 

Y dy Bo r r r' ’ 


*y. 




(A.l) 


(I - T ) 
' r r' 


and in the radiation relaxation region, see equations (9*3) , are: 


dT 


-sincp 


a 




To> ' 


3y « Va T a 3 ^/ 


3 p„ t (I- - * 


a a 


- j 

a a 


(A .2) 


The boundary conditions for equations (A.l) are: 

T 0 , I -*Oasy -♦ + » 

r } r 


and 


(A. 3) 


T 


0. 


J as y 

-r* * 'Y 


and the boundary conditions for equations (A .2) are: 


B i as y o - - «“ 


T -♦ 0_ , I VT , 

a B * a B "a 
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and 


(A .4) 


T -* 6 , I — J as y -* 0 

a a 7 a r ^a 


A variable transformation of d/dTi r = (l/T^) b/by^ and d/5T) a = 

(l/p T n ) b/by , respectively, to equations (A.l) and (A.2) makes 
a a a 

them more adaptable to numerical solution. After non- dime ns‘i onali - 
zation of the temperature by 0^, i .e . , t = T/Sg, and of the intensity 
by 0^, i.e., a) = 1/0^, equations (A.l) become: 


dut 

= _2 s 2 
dT 1 6 


5 

uT 


1 ’ 

K - \) 


(ui - T 4 ) 
r r ' 


(A. 5) 


subject to boundary conditions : 


\ 


-* + oo } a) - 0 as t -* 0 


and 


Tj 0 , u) -* J /0 as t 

'r J r r B r 


e r /9 B , 


and equations (A.2) beeoiie: 


(A.6) 


-i ,a (T a - i; 

dT a' 16 Cu) -ib ’ 
a a' 


dT} x , 

|a § 1 

dT £ 7 £7 * 

a tw - T ) 

a a' 


(A.7) 
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subject to boundary condition: 


T) - a. ai " l as t " i 
'a ’ a a 


and 


(A.8) 


V 0 ' “a J r /e B as T a “* 0 a /0 B = 0 r /0 B + 1 ' 


where $ = Bo sin<p/0g . 


Solutions to equation (A. 5) and (A. 7) subject to boundary condi- 
tions (A. 6) and (A. 8), respectively, are given in figures 5-8 for 
values of $ = 1.6, 4, 8 and l6 which correspond to values of Bo = 

0 . 2 , 0.5> 1.0 and 2.0 for sin<p = 1. 
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